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ABSTRACT 

We have observed the most distant (= 0.829) cluster of galaxies in the 
Einstein Extended Medium Sensitivity Survey (EMSS), with the ASCA and 
ROSAT satellites. We find an X-ray temperature of 12.3l'l2!2 keV this cluster, 
and the ROSAT map reveals significant substructure. The high temperature of 
MS1054-0321 is consistent with both its approximate velocity dispersion, based 
on the redshifts of 12 cluster members we have obtained at the Keck and the 
Canada-Prance-Hawaii telescopes, and with its weak lensing signature. The 
X-ray temperature of this cluster implies a virial mass ~ 7.4 x 10^^ M©, if 
the mean matter density in the universe equals the critical value (Qq — i), or 
larger ii Qq < 1. Finding such a hot, massive cluster in the EMSS is extremely 
improbable if clusters grew from Gaussian perturbations in an Qq = 1 universe. 
Combining the assumptions that f2o = 1 and that the intial perturbations were 
Gaussian with the observed X-ray temperature function at low redshift, we 
show that this probability of this cluster occurring in the volume sampled by 
the EMSS is less than a few times 10"^ Nor is MS1054-0321 the only hot 
cluster at high redshift; the only two other z > 0.5 EMSS clusters already 
observed with ASCA also have temperatures exceeding 8 keV. Assuming again 
that the initial perturbations were Gaussian and Q,o — 1, we find that each 
one is improbable at the < 10~^ level. These observations, along with the fact 
that these luminosities and temperatures of the high-^ clusters all agree with 
the low- 2; Lx — Tx relation, argue strongly that Qq < 1. Otherwise, the initial 
perturbations must be non-Gaussian, if these clusters' temperatures do indeed 
reflect their gravitational potentials. 

Subject headings: intergalactic medium - galaxies: clusters: individual 
(MS1054-0321) - X-rays: galaxies - dark matter - cosmology:observations - 
galaxies : interactions 

1. Introduction 

Clusters of galaxies occupy a special position among celestial self-gravitating structures, 
as they are the largest virialized objects in the universe. If all the cosmic structures we 
currently see grew hierarchically from small seed perturbations in the early universe, as 
many lines of evidence now suggest (e.g. Peebles 1993), then the largest clusters of galaxies 



-3- 



are also the most recent virialized structures to have appeared on the cosmological scene. 
The most massive of all clusters, which represent the most extreme high-density peaks on 
scales ~ 10^^ Mq in the initial field of density perturbations, correspond to several-cr events 
if the perturbations are Gaussian. 

Because massive clusters represent the high-cr tail of the distribution, their number 
density should rise dramatically as the overall amplitude of the perturbations grows. This 
characteristic of cluster formation renders massive cluster evolution particularly sensitive 
to the mean matter density of the universe, currently equal to Qq, in units of the critical 
density SHq/SttG. If ~ 1, perturbations should still grow like (1 + 2;)^^, and evolution in 
the number density of distant clusters with redshift should be rapid. If fig <^ 1, perturbation 
growth has virtually stopped, and massive cluster evolution should be much more mild. 

Evolution in the distribution of cluster masses sensitively indicates Qq, but measuring 
the virial masses of clusters (Myir) is difficult, even at low redshifts. A much more 
straightforward observational task is to measure how the X-ray temperatures {Tx) of 
clusters depend on redshift and then to use a model to relate Tx to Myi,. (e.g. Oukbir & 
Blanchard 1992). The Advanced Satellite for Cosmology and Astrophysics (ASCA) (Tanaka, 
Inoue, & Holt 1994) now makes such tests possible because it can measure the temperatures 
of the most luminous X-ray clusters at high redshifts (Yamashita 1994; Furuzawa et al. 
1994; Donahue 1996). Henry (1997) has recently published a temperature function for 10 
massive clusters at 2; ~ 0.3, chosen from the Einstein Extended Medium Sensitivity Survey 
(EMSS) (Gioia et al. 1990; Henry et al. 1992). His comparision of the ^ ^ 0.3 - 0.4 
temperature function with the ^ ~ temperature function suggests that Qq ~ 0.5 and rules 
out Qq = 1 with 99% confidence. 

In this paper we attempt to extract the maximum leverage out of this test for Qq by 
measuring the temperature of the highest-redshift cluster in the EMSS, MS1054-0321 at 
z — 0.83. We present both an ASCA spectrum and a ROSAT HRI image of this very 
luminous cluster (L^ = 5.6 x 10^ h~'^ ergs~^ 2-10 keV rest frame, — 1.1 x 10^^ ergs~^ 
bolometric). For qq — 0.1, — 7.7 x 10^^ ergs~^ 2-10 keV rest frame. Luppino & 
Kaiser (1997) provide a true color optical image and a weak lensing map. After presenting 
the observations in § 2, we examine how they constrain fio in § 3, finding that they strongly 
indicate < 1. Section 4 summarizes our results. The ASCA and ROSAT observations 
we present here are part of a larger program to measure the X-ray temperatures of three 
high-2;. X-ray luminous clusters in the EMSS, the other two being MS0451-03 at z = 0.53 
(Donahue 1996) and MS1137-I-66 ai z = 0.78 (Donahue 1998, in preparation). Throughout 
the paper, we parametrize Hq — lOO/i km/sec/Mpc and assume qo — 0.5, except where 
noted. 
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2. Observations and Data Reduction 

This section describes tlie data indicating that MS1054-0321 is an unusually massive 
cluster. We first discuss the ASCA spectra showing that MS 1054-0321 has a temperature in 
excess of 10 keV. Next we introduce the optical velocity dispersion and lensing data on this 
cluster, which agree with its high X-ray temperature. Then we show a ROSAT HRI image 
that reveals substructure in this cluster, suggesting that it may not be fully relaxed. 

2.1. ASCA Spectra 

ASCA executed a single long (~ 70, 000 s) exposure of the cluster MS1054-0321 
during May 23-25, 1995. This satellite carries four independent X-ray telescopes, two 
with Gas Imaging Spectrometers (GIS) and two with Solid-State Imaging Spectrometers 
(SIS), from which we obtained four independent datasets. The GIS detectors are imaging 
gas scintillation proportional counters with energy resolutions of 8% at 5.9 KeV and 
50' diameter fields of view. The two SIS detectors are X-ray sensitive CCDs with energy 
resolutions of 2% at 5.9 KeV and 22' by 22' square fields of view. These observations 
employed the SIS detectors in 2-CCD mode. However, we restricted the data analysis to 
the CCD containing the target and only used the second CCD to check the background 
estimates. 

To prepare the data for analysis, we extracted clean X-ray event lists using a magnetic 
cut-off rigidity threshold of 6 GeV c~^ and the recommended minimum elevation angles and 
bright Earth angles to reject background contamination (see Day et al. 1995). Events were 
extracted from within circular apertures of radii 3.25 arcmin (SISl), 2.53 arcmin (SISO), 
and 5.0 arcmin (GIS) to maximize the ratio of signal to background noise. We rejected 
~98% of cosmic ray events by using only SIS chip data grades of 0, 2, 3, 4, and we rejected 
hot and flickering pixels. Light curves for each instrument were visually inspected, and time 
intervals with high background or data dropouts were excluded manually. We rebinned 
the SIS data in the standard way to 512 spectral channels using Bright2Linear (see Day et 
al. 1995) with the lowest 13 channels flagged as bad, and we regrouped the GIS data so 
that no energy bin had fewer than 25 counts. Table |l| gives the resulting net count rates 
and effective exposure times. We do not expect the derived X-ray fluxes to be consistent 
between the detectors because the target was not centered and some of the extended flux 
may have missed the SIS detectors. 

Background estimates were taken from the regions of the detector surrounding the 
cluster detection, from summed deep background images supplied by the ASCA Guest 
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Observer Facility (GOF), and from the second CCD. The deep background spectra were 
extracted with the same spatial and background removal filters used to extract the source 
data. The intensities of the local background and the deep background were consistent with 
each other, but the uncertainties for the deep background estimate were larger (because the 
source extraction area was smaller). We therefore used the local background estimates for 
our analysis. 

We extracted SIS spectra binned in energy channels in two different ways: (1) binning 
by a factor 4 in energy, and (2) regrouping the data into energy bins containing a minimum 
of 16 counts. The first procedure provided slightly better constraints on iron abundance 
because it maintained some of the spectral resolution around the 6.7 keV Fe line, preserving 
the predicted contrast between the counts in the energy bin where the line should be 
and the continuum channels. The second procedure yielded slightly better contraints on 
temperature. After rebinning the data, we restricted our fitting of the SIS spectra to the 
0.5 — 5.5 keV range over which the signal-to-noise was adequate. The adequate range of the 
CIS spectra extended to 7 keV. 

To analyze the spectra, we used XSPEC (vlO.O) from the software package XANADU 
available from the ASCA GOF (Arnaud 1996). We fitted the spectral data from the four 
ASCA datasets and their respective response files (see Day et al. 1995). The individual 
SIS response matrices were generated with the tool sisrmg (1997 version), which takes 
into account temporal variations in the gain and removes the inconsistencies seen in data 
analyzed with the standard SIS response matrices. 

The standard model we used to fit the data was a Raymond-Smith thermal plasma 
with a temperature Tx, absorbed by cold Galactic gas with a characteristic column density 
of neutral hydrogen A^h- We did not detect a significant intrinsic absorption component, 
to a 90% upper limit of Nh < 1.4 x 10^^ cm~^. Thus, we fixed soft X-ray absorption 
at an HI value of 3.7 x lO^'' cm~^ (Gioia et al. 1990). Each SIS spectrum was fit with 
its own normalization; the normalizations of the CIS spectra were constrained to be the 
same. The parameters varied to provide a fit to 4 spectra were therefore the 3 independent 
normalizations (1 for each SIS spectrum, 1 for both CIS spectra), a single emission-weighted 
temperature, and the metallicity (iron abundance). Examples of the binned spectra and 
best-fit model convolved with each telescope's response matrix are found in Figure |^. The 
best-fit redshifted temperature to the binned-by-four spectrum was 12.4±|3keV with an 
upper limit on iron abundance of 0.22 solar, keeping the absorption fixed. If we used the 
SIS spectra regrouped to a minimum of 16 counts per bin, the best-fit temperature was 
virtually the same: 12.3±2 J keV with an upper limit on the iron abundance of < 0.25 solar. 
All uncertainties quoted are the 90% confidence levels for a 1-dimensional fit (Ax^ = 2.70), 
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Fig. 1. — We plot the X-ray spectra from each of the four ASCA telescopes for the cluster 
MS1054-0321. Each spectrum has over 1000 X-ray counts. (See Table 1.) The best-fit model 
is drawn with a sohd hne. The spectra we display were binned such that each energy bin 
contained a minimum of 16 (SIS) or 25 (GIS) counts. 
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and all fits have an acceptable reduced of ~ 0.8 — 0.9. Figure ^ shows the two-parameter 
contours for the cluster metallicity and Tx in keV. 

2.2. Optical Corroboration 

The available optical data on this cluster corroborate the high temperature measured 
by ASCA. Weak lensing observations of MS1054-0321 by Luppino & Kaiser (1997) show that 
its mass within 0.5 h^^ Mpc is (5 — 11) x 10^^ Mpc, depending on the redshift distribution 
of the lensed galaxies. In an isothermal potential, its mass distribution corresponds to an 
isotropic one- dimensional velocity dispersion ~ 1100 — 2200 km s~^, consistent with the 
value expected from the ASCA temperature: {kTx / fJ^rUpY^'^ = 1400 ± 170 km s~^ (90% 
confidence limits). 

In addition, we have obtained the redshifts of twelve cluster members with the CFHT 
Multi-Object Spectrograph (MOS) on April 13, 1994 and with the LRIS (Oke et al. 
1995) at Keck on January 4, 1995. These are listed in Table ^ Using Timothy Beers' 
program ROSTAT, we calculated their mean redshift and velocity dispersion, finding 
z = 0.829 ± 0.008 and a velocity dispersion of 1360 ± 450kms~^ (90% confidence interval). 
(The observed redshift dispersion has been corrected to first order for cosmological redshift 
by dividing by 1 + z). Several alternative methods (Beers, Flynn, & Gebhardt 1990) give 
the same answer within the uncertainties. While the small numbers of redshifts do not 
constrain the velocity dispersion very well, the value we find is consistent with both the 
lensing results and the X-ray temperature. 

Recent observations of other hot, high- 2; clusters further support the notion that their 
temperatures faithfully reflect their gravitational potentials. Figure ^ shows the relationship 
between the X-ray temperatures and velocity dispersions of three hot X-ray clusters at 
z > 0.5 from the EMSS and two cooler clusters that are companions of MS0016-I-06 at 
z = 0.54 (Hughes, Birkinshaw, & Huchra 1995; Hughes & Birkinshaw 1998). The X-ray 
temperatures for the three hot clusters are from this paper, Yamashita (1994), Furuzawa 
et al. (1994), and Donahue (1996), and their velocity dispersions are from this paper 
and Carlberg et al. (1996). The data for the two cooler clusters are from Hughes et 
al. (1995) and Hughes & Birkinshaw (1998). The hotter clusters agree with the relation 
kTx = fimpalj), while the cooler clusters might be slightly hot for their velocity dispersions, 
but nevertheless in agreement with the low-z cr-Tx relation of Edge & Stewart (1991a). 
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Fig. 2. — Wc plot the two-dimensional contours at 68.3%, 90% and 99% confidence levels 
(Ax^ = 2.30, 4.61 and 9.21) for the cluster iron abundance in units of the solar abundance 
and temperature in keV. These contours represent the confidence contours of simultaneous 
fits to the SIS spectra and to GIS spectra binned such that each bin contained a minimum 
of 16 and 25 counts respectively. 
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Fig. 3. — The relation between line-of-sight velocity dispersion ((Jid) and X-ray temperature 
{Tx) for several clusters at z > 0.5. The solid lines illustrates the relation kTx = fj,mpalj), 
and the other two lines show the ai^-Tx correlations derived from \ow-z clusters by Edge 
& Stewart (1991). The dotted line gives the correlation when Tx is the dependent variable, 
and the dashed line gives the correlation when (Tid is the dependent variable. Temperatures 
of high- 2; clusters at kTx ~ lOkeV appear to indicate the depths of their gravitational 
potentials with respectable accuracy. 
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2.3. X-Ray Imaging 

Our ROSAT data on this cluster suggest that we ought to interpret its temperature 
and velocity dispersion cautiously. We observed MS1054-0321 with the ROSAT HRI for 
over 121,590 seconds on December 6-18, 1996, yielding 1060 =b 90 net counts (including an 
earlier 12,136 second observation on May 5-6, 1996.) Figure § (Plate 1) shows a map of 
X-ray contours overlaid on a ground-based optical image of the cluster. For this figure, the 
X-ray emission was smoothed adaptively with Gaussians chosen such that the product of 
the number of counts and o"^ was roughly constant. The minimum contour corresponds to 1 
count per 8" x 8", which is ~ 3 — cr. The cluster emission is clearly extended over a scale of 
at least one arcminute. With the spatial resolution of the HRI, we have resolved the X-ray 
emission into two or three clumps and an extended component, clearly indicating that this 
cluster is not regular. 

King models centered on the cD galaxy in MS 1054-0321 tend to have high f3 values 
{P ~ 0.7 — 1.0) and large core radii (~ 1' or ~ 500h^Q kpc), although these values are 
not well constrained by the HRI data because of the cluster's irregular X-ray shape and 
dumpiness. A King model is rejected by the data, with reduced > 2 for both circular 
and elliptical King models. 

The most western peak is the most significant, containing ~ 144 ± 25 counts within 
24" above the smooth cluster component. The central peak is smaller, with ~ 90 ± 24 
counts and the easternmost peak is the least significant at 44 ± 14 counts above the smooth 
component. We made a basic statistical test for bimodality of the X-ray distribution of 
photons with the Lee statistic (Fitchett 1988; Lee 1979). We cahbrated the Lee statistic 
by producing multiple simulations of the null hypothesis in which the cluster emission is 
distributed in a King model (oc {r"^ + rl^re)'^^^^^'^) where the core radius is 56" and /5 = 2/3 
or 1. We found that the cluster is at least bimodal at greater than a 99.9% confidence level 
for either null hypothesis. 

The peaks in the X-ray emission do not coincide with any particular optical feature 
to within the 6 arcsec aspect uncertainty of the X-ray image. The cluster is elongated 
east-west, in the direction of the galaxy distribution, and a tail extends south of the cluster. 
Its X-ray emission roughly coincides with the dark matter contours as plotted by Luppino 
& Kaiser (1997), including the east-west and southerly extensions. The spatial resolution of 
the weak lensing map is not sufficient to confirm the sub-clumps seen in the X-ray map. The 
cluster X-ray emission also coincides with the shape of the Sunyaev-Zeldovich temperature 
decrement detected in centimeter interferometry experiments by L. Grego (PhD thesis, in 
preparation) . 
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The asymmetric appearance of the cluster in these high resolution images warns us 
that the relationship between its mass and its temperature might not be as simple as a 
spherical, hydrostatic, isothermal model would predict. However, recent hydrodynamic 
simulations of cluster formation beginning with cosmological initial conditions indicate that 
the systematic errors on cluster masses derived using a scaling relation between mass and 
temperature are relatively small, typically < 20% even when some substructure is present 
(Evrard, Metzler, & Navarro 1996). Higher resolution simulations of individual cluster 
mergers, whose initial conditions are not dictated by hierarchical structure formation, 
indicate that errors in the estimated mass can be somewhat larger 50%) after a recent 
merger, but these errors tend to yield underestimates of the cluster's mass (Roettiger, 
Burns & Loken 1996; Schindlcr 1996). The following discussion will therefore focus on 
the implications of the cluster temperature alone, assuming that the emission-weighted 
temperature is representative of the virial temperature. Mass errors at the 20% level do not 
affect the qualitative results of the ensuing discussion. 

3. Constraints on Temperature Evolution 

This cluster, MS1054-0321, is one of the hottest clusters known. If Qq = Ij structure 
formation simulations show that we can estimate its virial mass by assuming that the mean 
density within its virialized region is ~ 200 times the critical density at the cluster's redshift 
and that the cluster itself is isothermal (Evrard et al. 1996; Eke et al. 1996). The virial 
mass of MS1054-0321 would then be approximately My^, - [2kTx/i^mp{l + z)f/^{10HoG)-^ 
~ 7.4 X 10^'^ h"^ Mq within r2oo = l-5/i~^ Mpc. (We note that this relation is equivalent to 
assuming /3 = 1 in a standard X-ray mass relation for an isothermal gas.) Such a massive 
cluster at such a high redshift severely challenges models of hierarchical structure formation 
with ilo = 1 (Evrard 1989; Peebles, Daly, & Juszkiewicz 1989; Donahue 1993). This section 
elucidates how severe those challenges are. 

Because the virial masses of clusters are much harder to measure than their 
temperatures, we will restrict this analysis to focus on temperature evolution, relying on 
Tx to be a surrogate for Mvir. The analysis will therefore depend on the temperature- 
mass relation we assume. In an = 1 universe, this relation is relatively simple: 
Tx oc M^/'^(l + z). If i^o 1) this relationship breaks down below z ^ VIq^ — 1 a.s the 
development of structure stagnates. 
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3.1. Rarity of Hot Clusters at High- 2; 



In order to determine the rarity of a cluster as hot as MS1054-0321 in an f2o = 1 
universe seeded with Gaussian perturbations, we will assume that the Press-Schechter 
formula (Press & Schechter 1974) adequately describes the mass function of virialized 
objects, an assumption borne out by numerical simulations (e.g. Lacey & Cole 1994, Eke 
et al. 1996). In integral form, the comoving mass density of virialized objects with virial 
masses greater than M is then 



where po is the current mean matter density and Uc is the critical threshold at which 
perturbations virialize. This threshold is Uc = 6c{t)/a{M)D{t), where a{M) describes the 
rms amplitude of linear perturbations on mass scale M, D{t) describes their linear growth 
rate, and 6c{t) is the linear overdensity at which perturbations virialize. If Qq = 1, then 
z/c oc (1 + 2;) at a fixed mass scale M. 

This integral form of the Press-Schechter formula makes its theoretical underpinnings 
somewhat more apparent than the more familiar differential form does. The formalism 
assumes that the original linear density perturbations from which structure grew, smoothed 
over a length scale corresponding to mass M, follow a gaussian distribution with a dispersion 
proportional to a{M). As the perturbations grow, the most prominent of them collapse and 
virialize when their amplitudes exceed some threshold. The parameter i/c expresses how 
many standard deviations away from the mean amplitude a positive density perturbation 
must lie to have collapsed and virialized by time t. The total mass density in virialized 
perturbations of mass M or greater is then proportional to the integral over the tail of this 
Gaussian distribution from z/c to infinity, which yields the complementary error function 
in the above expression if o"(M) ^ as M — > 00. Because the typical amplitudes of 
perturbations decline with increasing M, the number density of virialized objects exceeding 
mass M is approximately M~^p{> M). At several standard deviations away from the 
mean, virialized objects with masses significantly exceeding M are much rarer than objects 
of mass M, so this approximation for the mean number density of virialized clusters is good 
to a few tens of percent. 

We will now proceed to determine the expected number density of clusters like 
MS1054-0321 at high redshift as follows: 

1. We conservatively assume that MS1054-0321 is hotter than 10 keV. 

2. In an Qo = 1 universe, Tx oc [1 -\- z) for clusters of a given mass, so the expected 
temperature at 2; = of a cluster as massive as MS1054-0321 is > 5.5 keV. 




(1) 
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3. Modifying the temperature-mass relation derived from simulations by Evrard et al. 
(1996) to reflect the mass within a region whose density is 200 times the critical value, 



This relation yields a virial mass > 7.4 x 10^^ h'^ M© for MS1054-0321. 

4. According to Henry (1997) the current number density of clusters hotter than 5.5 keV 
is < 3 X 10~^ Mpc~^. This number density decreases rapidly to higher temperatures 
and masses, so the mean virialized mass density in clusters hotter than 5.5 keV is 

< 2 X lO^h^ MqMpc-^. 

5. Invoking the integral Press-Schechter formula, we find that Uc{z = 0) > 3.4 for clusters 
with Tx > 5.5 keV. 

6. Devolving the growth of perturbations back to z = 0.83 implies that = 0.83) > 6.2 
for clusters with Tx > 10 keV. 

7. At z = 0.83, the mean virialized mass density of such clusters in comoving 
coordinates is thus < 3.1 x 10^ /i^ Mpc~^, and the corresponding number density 
is < 4.4 X 10-^3 Mpc"^ 

Next we will show that the detection of MS1054-0321 in the EMSS implies a space density 
~ 10~^ Mpc~^ for hke clusters, over 10^ times what we predict from the current cluster 
population and the assumption of Gaussian perturbations in an Qq = 1 universe. Before 
we go on, note that this test for is independent of the mass spectrum of perturbations 
because we are comparing clusters of similar masses (e.g. Oukbir & Blanchard 1997). 



The EMSS is particularly useful for studying cluster evolution because it has a 
quantifiable selection volume (see Henry et al. 1992; H92 hereafter). The usual procedure 
for determining the number density of clusters detected in such a survey is to compute T^iax 
the maximum volume within which a cluster of a given flux could be detected, for each 
cluster. The number density of high-redshift clusters in the survey is then the sum of the 
1/Vmax for all clusters in the redshift shell of interest. Here we consider the z = 0.5 — 0.9 
shell, which contains 4 EMSS clusters bright enough for ASCA to measure a temperature. 



we find 





3.2. 



Hot Clusters in the EMSS 
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Gioia & Luppino (1994) give the cluster X-ray fluxes within the detection cell for 
the EMSS clusters at z > 0.5. To determine V^^x for each cluster, we calculate the 
maximum redshift Zmax at which it could have been detected at each flux limit listed in H92, 
assuming that the cluster has a core radius of 0.25 Mpc and a /?— model surface brightness 
distribution with f5 = 2/3, while correcting for the fraction of the cluster emission that 
falls into the EMSS detection cell (H92; Donahue, Stocke & Gioia 1991). This correction 
factor is slightly smaller than two for z > 0.5 (H92). We then integrate over r'^dr with 
r = 2cHq^[1 — (1 + z)~^^'^] from the minimum redshift of the shell to either the maximum 
redshift of the shell or -Zmax? whichever is smaller. This computation gives the comoving 
volume per unit steradian associated with the cluster within the specified redshift shell in a 
flat universe. To determine Vmax we multiply this number by the solid angle corresponding 
to the appropriate flux limit. We compute V^ax assuming both (a) no /^-correction and (b) 
a i^'-correction with f,^ oc z/~°'^. Varying any of these assumptions within reasonable ranges 
changes the results by ^ 10%. 

The Knax associated with MS1054-0321 in a flat universe is 6.9 x 10^ h'^ Mpc^ with 
no i^'-correction and 6.5 x 10'' h'^ Mpc'^ with the i^-correction. These volumes imply a 
comoving space density ~ 1.5 x 10~^h^Mpc~^ for clusters like MS1054-0321 in the redshift 
shell z = 0.5 — 0.9. Over four orders of magnitude separate this density from the expected 
value ~ 4 X 10~^^ /i^ Mpc~^, casting severe doubts on either the assumption of = 1 or 
the assumption of Gaussian initial perturbations. 

While this single hot high-z cluster argues strongly against r2o = 1) it does not 
constitute an airtight case on its own. One could argue that MS 1054-0321 is anomalous 
in some way, but there are two other high-2 clusters in the EMSS that are nearly as 
hot: MS0016+16 with Tx = 8.4 keV at z = 0.54 (Yamashita 1994; Furuzawa et al. 
1994) and MS0451-03 with Tx = 10.4 ± 1.2 keV at z = 0.54 (Donahue 1996). Figure | 
demonstrates the consistency between their X-ray temperatures and optical velocity 
dispersions. (MSI 137-1-66, the fourth high redshift EMSS cluster, was observed by ASCA 
in 1997, and will be discussed in a separate paper.) 

For comparision, we can compute the expected number density for these clusters 
based on the same assumptions as before. Given Tx > 8 keV at z = 0.5, we find 
Mvir > 6.9 X 10^^ Mq, corresponding to Tx > 5.3 keV at z = 0. The current mean 
virialized mass density in such clusters is similar to the case of MS1054-0321, implying 
Pc > 3.4 at z = and z/^ > 5.2 at z = 0.54. For fig = 1 and Gaussian perturbations, the 
expected number density of > 8 keV clusters at z = 0.54 is then < 8.6 x 10~^^ /i^Mpc"^. 
Each of the two hot EMSS clusters at z = 0.54 would therefore be rarer than one in a 
hundred, with a joint improbability of order a few times 10~^, somewhat more probable 
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than the chance occurrence of MS1054-0321. 

In summary, given the three hot high redshift clusters in the EMSS, the imphed 
number density of clusters with Tx > 8keV at 2; = 0.5 - 0.9 is ~ 1.6 x 10"^ Mpc~^. The 
probability of all three clusters appearing in the volume sampled by the EMSS in an = 1 
Universe is ~ 10"^ if the initial density perturbations have a Gaussian distribution. 

Universes with ~ 0.3 have no problem accomodating such hot clusters at these 
redshifts. Eke et al. (1996; see also Viana & Liddle 1996) show that the cluster temperature 
function evolves very little in an = 0.3 universe with no cosmological constant and only 
modestly in a flat universe with the same Q,q. The number density of > 8 keV clusters from 

Henry (1997) is ~ 5 x 10~* /i^ Mpc~^, quite close to the value we find at z ~ 0.5. The 
statistics of hot clusters at 2; > 0.5 therefore add substantial weight to the growing body of 
evidence emerging from cluster evolution that VLq ^ 0.3 — 0.5 (e.g., Carlberg et al. 1997; 
Henry 1997; Bahcall, Fan, & Cen 1997). 



3.3. High-z Lx-Tx Relation 

Cluster evolution from z ~ 0.5 to the present should be modest at best in a universe 
with low VLq. Here we show that the relationship between the X-ray luminosity {Lx) and the 
X-ray temperature of MS 1054-0321 is similar to the relationship derived from lower- redshift 
clusters. Edge & Stewart (1991b) found that Lx ~ (8.8 x 10^^ /i'^ erg s-^){Tx/l keV)2-79±o.o5 
in the 2—10 keV band for low-redshift clusters, while Mushotzky & Scharf (1997) demonstrate 
that clusters observed at 2; > 0.14 by ASCA, including MS0451-03 (Donahue 1996) and 
MS0016-I-16 (Yamashita 1994; Furuzawa et al. 1994), do not deviate appreciably from 
this relation. Our observations indicate that Lx{2 — lOkeV) ^ 7.7 x 10^^ erg s~^ 
for MS1054-0321 in its own rest frame {qq — 0.1), implying an expected temperature of 
11.3 keV, on the lower edge of our one-dimensional 90% confidence interval. 

The lack of evidence for Lx-Tx evolution in our observation of MS1054-0321 is 
consistent with the other evidence for < 1, but this test for is not as powerful as 
the temperature function test, as discussed in Section 3.2. Some cluster evolution models 
predict little evolution in the Lx-Tx relationship in an f2o = 1 universe, as long as the 
minimum entropy of the X-ray gas remains the same for all clusters regardless of redshift 
(Kaiser 1991, Evrard & Henry 1991, Bower 1997). The lack of evolution in Lx-Tx would 
then require some sort of non-gravitational heating of the X-ray gas. 
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4. Summciry 

The ASCA observations we have presented here imply that the EMSS cluster 
MS1054-0321 at z = 0.83 has an X-ray temperature of I2.3I2J keV. Supporting optical 
redshifts of 12 member galaxies (this paper) and the weak lensing map (Luppino & Kaiser 
1997) provide a velocity dispersion and a lensing mass that agree with the mass inferred 
from the high X-ray temperature. Our ROSAT image exhibits significant substructure, but 
numerical modelling of cluster evolution has shown that the presence of substructure does 
not preclude using a cluster's X-ray temperature to estimate its virial mass (e.g. Evrard et 
al. 1996). 

This high a temperature in so distant a cluster is powerful evidence that Qq < 1- Using 
the integral Press-Schechter relation and the current temperature function of clusters, 

we calculate the expected number density of such hot clusters at z = 0.83, under the 
assumptions of Gaussian perturbations and Qq = 1. We find that the product of the 
expected number density and the effective volume of the EMSS in the z = 0.5 — 0.9 redshift 
bin is less than a few times 10~^. Similar calculations for two other hot EMSS clusters at 
z ~ 0.54 give probabilities less than a few times 10~^ for each one. 

For Qq ^ 0.3, the temperature function of hot clusters is expected to evolve very little 
to redshifts > 0.5 (Eke et al. 1996), so that the hot EMSS clusters are quite consistent 
with a low-^o universe. We also find that MS1054-0321 luminosity and temperature are 
consistent with the low- 2; Lx-Tx relation, a result also consistent with Qq < 1. These 
results strongly support the mounting body of evidence (e.g. Carlberg et al. 1997) showing 
that the mean density of matter in the universe is fig ~ 0.3. 
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Table 1: ASCA Observing Information 
Detector Aperture Radius Count Rates Useful Exposure Time 



Arcniin s " s 

5150 2^53 0.019 62,019 

5151 3.25 0.013 58,176 
GIS2 5.00 0.012 65,534 
GIS3 5.00 0.016 65,920 



Table 2: Member Galaxy Redshifts 

ID Redshift Offset E/W (+/-") Offset N/S (+/-") 



(cD) 


0.8309 ±0.0008 


0.0 


0.0 


2 


0.8342 ±0.003 


+23 


-6 


3 


0.8127 ±0.0003 


+56 


-5.5 


4 


0.8213 ±0.0007 


+37 


+30 


5 


0.830 : ±0.008 


-14 


-20 


6 


0.8209 ±0.001 


-29 


-13 


7 


0.8286 ±0.001 


-31.5 


-11.5 


8 


0.8353 ± 0.0006 


-38 


-8 


9 


0.8332 ±0.001 


-43 


-6 


10 


0.8143 ±0.002 


-38 


+0.5 


11 


0.8378 ±0.003 


-80 


-44 


12 


0.8319 ±0.002 


-97.5 


-38.5 
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Fig. 4.— Plate 1. The central 3.75' x 3.75' of the cluster MS1054-0321. The optical image is 
a 14,400 second I-band image taken with the University of Hawaii 88-inch (Gioia & Luppino 
1994). The ROSAT HRI image, a sum of all available HRI data for a total livetime of 121,590 
seconds, was rebinned into 8" x 8" pixels and adaptively smoothed with three Gaussians with 
sigmas of 16, 14.3, and 12.6 arcseconds corresponding to counts in a pixel of < 10, 10 — 13.5 
and > 13.5. The sigmas were chosen so that the products of the number of counts and 
were roughly constant. The background was 9.4164 counts per pixel. The net contour levels 
are 1, 2.8, 4.6, 6.4, 8.2, and 10 counts per pixel. These net count rates correspond to X- 
ray surface brightnesses of 1.4, 3.9, 11.6, and 14.0 x 10~^^ erg cm~^ s~^ arcmin"^ assuming 
that 1 HRI count ~ 2.8 x 10~^^ erg cm~^. The galaxies are marked with identifications 
corresponding to Table |^. The central galaxy is at RA(2000) = 10 56 59.9 and Dec(2000)=- 
03 37 37.3, as measured from our HST image of this cluster (Donahue et al., in preparation). 




06^ 04^ 02^ 58^ 56^ 54^ 

10"^ 57 "^00^ 



